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Poa Wee 


Experimental investigations have been conducted to 
examine the general problem of convective heat transfer ina 
rectangular cross section channel with both straight and 
curved passages. The experimental configuration was modeled 
as infinite parallel plates with one wall at a constant heat 
flux and the opposite wall adiabatic. The air flow through 
the channel was heated through the Joulean heating of a wall 
composed of Temsheet. Liquid crystals applied to the surface 
of the Temsheet were used to confirm the presence of the 
Taylor-Gortler vortices. Average Nusselt number measurements 
were used to investigate their effect on heat transfer from 
the curved test section. Experiments were conducted for 
steady state laminar flow that was hydrodynamically developed 
and thermally developing. It was confirmed that the Taylor- 
Gortler vortices enhanced the heat transfer process in the 
curved section. The experimental data was compared with 


several analytic studies with good results. 
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A. DESCRIPTION OF TAYLOR-GORTLER VORTICES 

A considerable amount of data and indirect evidence, [1, 
ee 3 | among others, has shown that the fully developed laminar 
flow of viscous fluids along a concave wall does not remain 
two dimensional. Instead, the flow forms spiral vortices in 
counter-rotating pairs which frequently possess a regularly 
spaced cellular structure. This phenomenon, known as 
Taylor-Gortler vortices, occurs as a result of flow insta- 
bility induced by variations in the centrifugal forces acting 
on fluid particles. The Taylor-Gortler vortices are laminar 
in nature with their axes aligned in the eon of the 
principal fluid flow and with secondary flow velocities in 
both perpendicular directions. An example of this flow 
Pacem 1S Shown in Figure 1. 

In a channel that is curved in the streamwise direction, 
such as that used in this study, the fluid particles near 
the center are subjected to higher centrifugal forces than 
those experienced by the slower moving fluid particles near 
the boundary of the concave wall. Consequently, the tendency 
is for the fluid particles in the center of the channel to 
move outward to the concave wall. Under this influence, 
particles near the wall are displaced initially in the 


Spanwise direction and subsequently continue radially inward, 
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Peplacing the central fluid particles. Once in the center 
Ofethe channel the fluid particle will come under the in- 
fluence of the higher free stream velocity and resultant 
centrifugal forces. This action forces movement back toward 
the concave wall causing the entire cyclic rotation process 
me propagate. 

Although the Taylor-Gortler vortices are laminar, it is 
quite possible that they play an important role in the transi- 
tion from the laminar to the turbulent flow regime. [2] 
Additionally, the presence of the vortices have been thought 
to account for the increase in heat transfer from curved 
walls due to mass transfer. [4] 

The Taylor-Gortler vortices also have similarities to 
other vortex flow patterns such as the longitudinal vortex 
rolls developed in the laminar forced convection heating 
of fluid layers between parallel plates. [5] In another 
field it has further been hypothesized that the cross-hatching 
phenomenon observed in re-entry vehicles can be explained 
by the presence of centrifugally induced secondary motion 
characteristic of the Taylor-G6rtler vortices. [6] 

It is apparent that the understanding of these vortices 
and their effect on heat transfer and fluid flow character- 
istics could significantly effect many areas of fluid 
mechanics. Recently it has been proposed es a6 g | that 
concave wall curvature can have an advantageous effect on 


film cooling over turbine blades. 
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ie spp eatieons Of =lurid Mechanics and heat transfer 
[mmvolve curved surfaces or curved channels and the under- 
standing of the fluid mechanics in this geometry will greatly 


enhance engineering designs and capabilities. 


Epo RLEF HISTORY 

ijstability Of an inviscid fluid in a curved flow was 
first considered by Lord Raleigh in 1916. [10] He concluded 
meecche Stability criterion for an inviscid fluid had as a 
necessary condition that the circulation always increased 
Wee the radius. In essence, this meant that the product 
of the local circumferential velocity and the corresponding 
local radius of curvature cannot decrease with increasing 
radial distance. 

iieoeomcem!. Taylor expanded the scope of Peete 
meet y criterion to include viscous fluids. [1, 11] His 
extensive analytic and experimental investigation focused on 
the flow of viscous fluids between a rotating inner cylin- 
der and a stationary outer cylinder. Taylor ascertained that 
such couette flows become unstable when a characteristic, 
non-dimensional parameter exceeded a critical value of 4l. 
Referred to as the Taylor number, this parameter was defined 
as: 
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In this relationship d represents the spacing between the 
concentric cylinders, assumed small when compared to the 
radius of curvature, R. Re is the Reynolds number based on 
U, the circumferential velocity of the inner cylinder and d, 
the annulus width. When flow conditions were such that the 
Taylor number exceeded this critical value, the secondary 
flow observed in the gap between the stationary outer cylin- 
der and the concentrically rotating inner cylinder was in the 
form of longitudinal vortices whose axes were aligned in the 
circumferential direction. These vortices were given the name 
itamalor vortices. 

A similar form of instability occurs when a viscous 
fluid flows in a curved channel due to a pressure gradient 
acting along the concave wall. W. R. Dean m2] first addressed 
thas problem analytically in 1928. Specifically, he consi- 
dered flow in a channel formed by two concentric cylinders with 
the restriction that the cylinder gap distance be small com- 
pared to the radius of the inner cylinder. Dean concluded 
that the flow instability would initiate, and the secondary 
Micnewertices (Similar to the Taylor vortices) would form, 
when the Dean number was greater than 36. The Dean number 


Wee then defined as: 


wa 
re J 


De 
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In this case, the Reynolds number was based on the mean velo- 
city of the undisturbed flow and d represented the channel 
half width. The analytical work by Dean was later verified 
by W. H. Reid [13]J using an approximate numerical solution. 

ils Oesteecorticr [2] studied the influence of small 
disturbances on the stability of laminar boundary layer 
profiles. Gortler found that these disturbances closely 
resembled those described by G. I. Taylor which led to fluid 
Meee mstabllity in the form of vortices. In his numerical 
calculations, Gortler concluded that the amplified distur- 
bances were produced only on the concave walls; nonetheless, 
the overall flow regime seemed to remain laminar. The appro- 
ximate analytic results that Gortler had developed were 
initially verified with an exact solution by G. Hammerlin as 
reported by H. Schlicting [14], and later substantiated by 
BEMCXEeNSiVe Numerical analysis by A. M. O. Smith [3]. M. 
Pewtenersset ale, [15] and S. Winoto, et al., £16], among 
others, have been able to verify these numerical predictions 
by using hot wire anemometry, laser Doppler systems and flow 
visualization techniques to describe the motion of the secon- 
dary flow. A non-linear analysis of Gortler vortices was 
eemauercd by Y. Aihara {17] in 1976. 

Growing interest in the various effects of the secondary 
Flow associated with Taylor-Gortler vortices led to more 
specific studies concerning the vortice influence on heat 


Mamita mmine both turbulent and laminar flow. In 1955, 
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F. Kreith [4] reported that the heat transfer from a heated 
concave wall was considerably higher than that from a convex 
Toeieor Ladentical curvature under similar conditions of tur- 
bulent flow. 

In the special cases of very high or very low Prandtl 
mumperc, Lb. Persen [18] related the increased heat transfer 
from curved walls to the presence of Taylor-Gortler vortices. 
He determined that their overall effect was to enhance heat 
transfer through the boundary layer. 

In 1970, McCormack, et al., [19] published the results 
of their initial experimental work regarding the effect of 
secondary flow on heat transfer. From the data obtained they 
concluded that a complete theoretical explanation of the heat 
transfer effects required that all non-linear terms be re- 
tained in the flow equations. R. Kahawita and R. Meroney [20] 
considered this hypothesis and further reported that both 
the high order terms and the normal velocity components of 
the main flow were essential to any calculations at small 
wave numbers. 

i@uecGione Mort and Y. Uchida £5) considered the fully 
developed forced convection heat transfer between horizontal 
plates. In the results of their investigation they reported 
that as the temperature difference between the plates was 
increased a critical value developed beyond which longitudinal 


vortex rolls, similar to the Taylor-Gortler vortices, would 
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Remmewith their axes in the streamwise direction. M. Akiyama, 
et al., [21] confirmed the presence of these longitudinal 
vortices at a critical Rayleigh number of 1708 in cases where 
the plate is heated from below. 

Although heat transfer in curved duct flow has, and will 
continue to be investigated experimentally, considerable 
emphasis has been directed recently towards heat transfer 
predictions using numerical approximation for flows in curved 
duct geometries. In 1970, K. Cheng and M. Akiyama [22] 
developed a numerical solution for forced convection heat 
transfer in laminar curved rectangular channels by using a 
point successive over-relaxation method. For similar flow 
and boundary conditions, G. Yee and J. A. C. Humphrey [23] 
used the fully elliptic forms of the transport equations to 
compose a numerical solution. However, both of these numer- 
ical solutions were limited to small aspect ratios of essen- 
tially square cross section. 

In contrast to the large volume of published numerical 
approximations to this heat transfer problem, the number of 
experimental and analytical investigations have been very 
limited. Results for a fully developed, constant wall heat 
flux, curved, square cross section channel were reported by 
Yo. Mora. ec al. [24], and results for larger aspect ratios 


were reported by M. Durao. [25] 
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i Le eon THE SlUuby 


In view of the possible engineering design implications 
of the Taylor-Gortler vortex phenomenon, the primary purpose 
of this study was to investigate the effects of the secondary 
flow on the heat transfer rates in a curved channel of rec- 
Gemewrar cross section. 

It was expected that the presence of the Taylor-Gortler 
vortices would result in the secondary flow velocity com- 
ponents transporting the heated fluid particle from the cur- 
ved wall to the center of the flow channel. This fluid 
motion would displace the cooler fluid particle in the center 
of the channel causing it to move toward the heated curved 
wall. It was anticipated that this cyclic motion would 
measurably improve the heat transfer process for all fluid 
velocities in the laminar flow regime. Secondary investiga- 
tions were conducted to determine the additional influence 
of buoyant effects on the development of heat transfer in 
the curved channel. It was expected that the heat transfer 
rates would be greater when the buoyant forces were in the 
direction of the main flow. 

The investigation was conducted using a rectangular cross 
section channel that contained both a straight and a curved 


test section. The results obtained in the straight 
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section served as a baseline and as such were compared to 
meemcata Obtained im the curved section. The complete 
process was accomplished for two distinct channel orienta- 
tions: (1) buoyant forces aligned with, and (2) buoyant 
forces opposed to, the main direction of the fluid flow. 
Straight section data were compared with analytic solu- 

tions provided by P. McCuen, et al., [26], P. M. Worsde- 
pemmras [27] and R. K. Shah and A. L. London. [28] The 
buoyant effects were qualitatively compared with the data 
contained in a numerical computation by R. Chilukuri and 


J. A. C. Humphrey. [29] 
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Pee DESCRIPTION OF THE APPARATUS 

A rectangular cross section channel built as described 
in reference 25 was used to achieve the objectives of this 
study. The channel was manufactured from two 0.635 cm thick 
sheets of Plexiglas separated by 0.635 cm thick spacers whose 
inside lengths served as the sides of the channel. 

As shown in Figure 2, the channel was composed of a 
Peetont portaenm 122.0 cm long, an intermediate curved 
section which formed a 180 degree are and a final short 
straight section. The radius of curvature of the interior 
concave wall of the curved portion was 30.5 cm. The channel 
was 0.635 cm high and 25.4 cm wide with an aspect ratio of 
Momadd a cross sectional area of 16.13 square centimeters. 

A eross sectional view of the channel is shown in Figure 3. 

The working fluid was air at room temperature which 
entered the flow channel through an entrance bell constructed 
of Plexiglas. The entrance bell was designed and manufactured 
in accordance with ASME nozzle standards with an elliptic 
curve based on a ten inch major and one inch minor axis. A 
permeable cloth covered the entrance nozzle to reduce the 
introduction of dust and other foreign matter into the flow 


channel. The flow at the end of the channel exited through 
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an aluminum exhaust bell connected by flexible tubing to a 
Fisher and Porter Company variable area flow meter, model 
number 10A3565A. The rotameter had a 100% full scale rating 
ol cubge meters of air per minute (11.1 standard cubic 
feemeot air per minute) at one atmosphere and 70° F. 

The flow of air was drawn through the channel, the flexi- 
ble tubing and the rotameter by a single electrically driven, 
Cadillac centrifugal blower, model G-12. The blower speed 
was controlled by varying the motor voltage with a General 
Radio Company Variac Autotransformer type WILOMT3. 

For the purpose of obtaining experimental heat transfer 
fees two test sections were Constructed; (1) a straight 
fesmesection, 29.2 ecm im length, located in the straight 
portion of the channel downstream of the hydrodynamic entrance 
region, and (2) a curved test section 28.3 cm in length, sub- 
tending a circular arc of 53.1° in the lower quadrant of the 
curved portion of the channel. The area of the straight 
test section was 741.9 square centimeters and that of the 
curved test section was 717.7 square centimeters. 

In each test section the outer wall was detached and 
modified as shown in Figures 4 and 5. Temsheet, which is a 
carbon impregnated porous paper with the property of uniform 
electrical resistivity, was glued to the entire interior 
surface of the outer wall in each section. The flow of air 
could then be heated through the Joulean heating of the Tem- 
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was slightly temperature dependent, a precision resistor with 
pe lectrical resistance of 2.0236 ohms was connected in 
series with the Temsheet to facilitate calculating the 
instantaneous power supplied. A LAMBDA Regulated D.C. Power 
Supply model LK345A FM supplied electrical current to heat 
the Temsheet to desired levels. 

The variables to be measured in the investigation were: 
the flow inlet temperature at the entrance to each test 
section (Tey? > the flow outlet temperature at the exit of 


each test section (T ), the wall temperature of the heated 


OUT 
plate in each test section ae the wall temperature oppo- 
Site the heated plate in each test section (T 3)> the flow 


exit temperature at the end of the channel (T ), the temper- 


eee 
atures between the two insulation layers for each test section 


(T ) 


Ins?? and the voltages across the precision resistor (Vi5e 
and the heater (Vi). Copper-Constantan thermocouples were 


used to measure all the above temperatures. 


midiid Crystals 


Cholesteric Liquid Crystals were employed to provide 
Mieiehecdemyc Meaisirements Of the surface temperature distri- 
bution on the Temsheet. The liquid crystals exhibit drama- 
tic changes in color for a small temperature differential. 
The liquid crystals would allow a reasonable qualitative 
comparison to be made between the heat transfer processes on 


each test section by observation of the color pattern 
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displayed for various temperature gradients and flow velo- 
@9ties. opecifically, the occurrence of multi-color stripes 
aligned in the flow direction on the curved wall would indi- 
cate regions of high and low temperatures and confirm the 
presence of Taylor-Gortler vortices. 

To select the appropriate liquid crystal system for 
application to the surface of the Temsheet, consistent 
values of flow inlet temperature, Reynolds number, Prandtl 
number and Nusselt number were assumed and an approximate 


energy balance was developed: 








Q = mC oair‘Tour ~ brn? = 2 Ap, AT 
Or 
aol z Re Pr (A) 
Tour ~ Tin Nu Apr) 
where 
ee. ON Our 
wo Z 


A desired ten degree Celsius difference between the 
inlet and outlet temperatures required an average wall 
temperature of forty-five degrees Celsius. Taking into 
consideration the possible losses in the system, the required 
power to achieve that average wall temperature was estimated 


to be a maximum of 40 Watts. 
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Eponecompletion Cf a series of preliminary experiments it 
Peieowestablished that a combination of liquid crystals would 
be required to adequately delineate the average wall tempera- 
ture and the vortex phenomenon. The test sections were 
therefore sprayed with a mixture of three NCR Microencapsu- 
lated Liquid Crystals (R-41, R-45 and R-49) prepared in 
accordance with reference 30. Using an artist's air brush, 
approximately twenty coats of the liquid crystal mixture were 
applied directly to the Temsheet in order to ensure a sharp 
color contrast and definition within the response curve of 
each liquid crystal. It was hoped that the resultant color 
patterns could be recorded photographically. 

A total of sixty-four copper-constantan, glass insulated, 
30 gauge thermocouples were constructed and placed in specific 
locations to measure the required temperatures. All ther- 
mocouples were permanently wired to a Hewlett-Packard Data 
Acquisition System, Model 2010C, which provided an automatic 
printed record of the instantaneous data for each thermo- 
couple position. The full set of thermocouples were connected 
in series to a Kaye Instruments Ice Point Reference, Model 
K-140-4. In each test section five thermocouples were inser- 
ted between the first and second layers of the insulation 
with an additional five inserted between the second and third 
layers. Each set of five thermocouples were connected in 


parallel to read an averaged temperature value (Trig? between 
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the respective insulation layers. A set of four thermo- 
couples, also connected in parallel, were inserted into the 
flow channel at each of the following locations: the en- 
trance to the channel, the exit of the straight test section, 
and the exit of the curved test section. At each location 
the thermocouples were distributed in a spanwise interval of 
eeoecm and at a d@épth of 0.3 cm to allow the thermocouples 

to read a bulk temperature for the air flow. In each test 
section eight thermocouples were placed in direct contact 
with the Temsheet through very small diameter holes drilled 
in the Plexiglas above the Temsheet and the thermocouple 
beads were electrically insulated from the Temsheet with 
ENMAR Heat Resisting Glyceryl Phthalate. An additional 

eight thermocouple beads were embedded in the Plexiglas 
channel opposite those positioned in the Temsheet. Figure 

6 provides a detailed sketch of this thermocouple arrange- 
ment. The thermocouples in contact with the Temsheet pro- 
vided an average outer wall temperature i ae and confirmed 
temperature distributions indicated by the liquid crystal 
patterns. Those in the Plexiglas channel wall provided 

an average inner wall temperature (T and assisted in 
estimating the losses in the system due to radiation. All 
the thermocouples were calibrated using a ROSEMOUNT Communta- 
ting Bridge, model 920A, and a ROSEMOUNT Constant Temperature 


Bathn Model 913A. 
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The thermal insulation consisted of Armstrong ARMAFLEX 
Bees Iisttarion which is a flexible foam plastic material. 
Three layers of insulation were used, each one with a 
thickness of 0.635 cm and covering an area slightly larger 
than each test section. The insulation was fixed in place 
with the use of an Adhesive Heat Resistant Ventilation Duct 
Tape as shown in Figure 7. The entire length of the channel 
was sealed with GENERAL ELECTRIC Silicone Rubber sealant 
caulk. 

Two aluminum electrodes, 0.318 cm thick, were inserted 
between the Temsheet and the Plexiglas as shown in Figure 3. 
Fach plate consisting of the Temsheet, the liquid crystals, 
the thermocouples and the aluminum electrodes was affixed in 
the channel using Teflon screws to provide electrical insula- 


teal T) « 


Beene > RIMENTARSEROCEDURES 

Experiments were conducted for each test section following 
the same set of procedures. The data obtained in the straight 
test section provided baseline information for later curved 
section runs, either aligned with or opposed to the buoyancy 
forces. The volumetric flow rate was varied incrementally 
from 0.063 to 0.220 cubic meters of air per minute to corres- 
pond to values of Reynolds numbers in the laminar flow regime. 
Data was recorded and calculations were performed for each 


FLOW Pateminerement . 
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The appropriate electrical power required to heat the 
Temsheet to a temperature level that would initiate the 
liquid crystal response was determined through preliminary 
experiments. 

Estimates of the time required for the experiment to 
reach steady state, decisions on what data would be recorded, 
and frequencies at which that data would be obtained were also 
based on preliminary runs. 

The determination as to what constituted a steady-state 
condition required careful consideration in light of the in- 
fluence that the multiple input variables could have on the 
results. The final criteria for steady-state was based on 
three input variables; the Flow rate percentage, the outside 
wall temperature (Temsheet), and the heater voltage. When 
these input conditions had a variation of less than two 
percent over a ten minute interval, the experiment was consi- 
dered to have reached steady-state. Under these guidelines 
the assumptions of steady flow and constant heat flux from 
the upper wall were valid. 

The time required to reach steady-state was found to be 
between two and two and one-half hours. Accordingly, all 
runs were two and one-half hours in duration. Data was 
recorded at one hour intervals for the first two hours and at 


ten minute intervals thereafter. 
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The instantaneous values for all the temperatures were 
peo matcically printed by the digital recorder in milivolts. 
A standard thermocouple table was used in conjunction with 
previously obtained thermocouple calibration data to convert 
the milivolts to degrees Celsius. The volumetric flow rate 
was read directly from the rotameter. The precision resis- 
tor voltage (Vip? and the heater voltage CV.) were recorded 
in volts and the instantaneous voltage supplied so was 


calculated using the following relationship: 


where Rop was 2.0236 ohms, the measured electrical resistance 
Smee precision resistor in the circuit. 

The steady state response of the liquid crystals was 
recorded photographically using a NIKOMAT (FIN) 35mm single 
reflex camera with a 85 - 205 mm zoom lens. The film was 
25 ASA Kodachrome with the shutter speed and lighting con- 
trolled manually. 

All data reduction was accomplished by desk-top computer 


programs and stored routines. 
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Cece oenN eA iON OF DATA 


A. ANALYSIS 

The flow examined in this study was laminar and steady, 
thus in a general sense, plane Poiseuille flow. The Temsheet 
which possessed the property of uniform electrical resistivity 
approximated a constant axial and peripheral wall heat flux. 
Temperature measurements on the Plexiglas’ wall opposite the 
heated surface indicated negligible heat loss through that 
wall, 1.e., adiabatic. With an aspect ratio of 40 the chan- 
nel was approximated by infinite parallel plates. The 
straight portion of the flow channel, upstream from the 
straight test section, was of sufficient length to ensure 
that the flow was hydrodynamically developed for all flow 
velocities in the laminar range. For the purposes of this 
study, the flow was considered laminar when the value of the 
Reynolds number was less than 2000. Based on the above 
assumptions, the experimental configuration was then modeled 
as laminar forced convection between parallel plates with 
hydrodynamically developed and thermally developing flow. 
The boundary conditions are fundamental boundary conditions 
of the secondkind with one wall at a specified axial and 


peripheral heat flux and the other wall insulated. [28] 


oy 





The appropriate variables, corresponding equations and solu- 

tion procedures are contained in the following paragraphs. 
The average heat transfer coefficient from the heated 

wall to the flow of air was defined by a form of Newton's 


ia OL cooling: 


Qair a ou a 


iieciis expression, — was the heat convected to the air, 
h was the average heat transfer coefficient, A,, was the 
surface area of the heated Temsheet, and AT was the tempera- 
ture difference between the average wall temperature Cae 
and the bulk temperature of the fluid (T,)-. 


The actual heat convected to the air ee was then 


independently calculated with the equation: 


oe eo mC (T 


eee) 
air pair 


OUT “IN 


where eae was the specific heat of air at constant pressure 
and m was the mass flow rate of the air. The mass flow rate 
was the product of the volumetric flow rate (q) and the local 
Flow density in the rotameter (9). The volumetric flow rate 
was corrected to standard air conditions for the rotameter 
and the local flow density was calculated assuming perfect 
gas behavior at constant pressure and using the flow exit 


De 


temperature (Toyter 
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The average Nusselt number was determined from the aver- 


age heat transfer coefficient using the following expression: 





Nu 


where d represented the height of the channel and ae was 
the thermal conductivity of the air evaluated at the bulk 
temperature (T2). 


The Reynolds number was calculated for all the experimen- 


maleruns and was defined as: 


_ pud _ md 
Re = sae uA, 


where d was the height of the channel, uw was the dynamic 
viscosity of the air at the flow exit temperature and A, 
was the cross sectional area of the channel. 


For all experiments in the curved test section, the Dean 


mnumber was calculated as: 


De = Re v at 


il 


In this equation d was again the height of the channel and 
Reewas the radius of curvature of the inner convex wall. 


To verify that Q which was indirectly measured by the 


ain: 
temperature rise across the channel, was actually the con- 
vection heat transfer from the heated Temsheet surface, an 


Poin neem radtarion heat transter was made for 6ach run. 


The average temperature of the lower plate was measured by 
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6ight thermocouples and the radiated heat transfer (Q\. was 
mmenlated with the Stefan-Boltzman law of radiation: 
Ci ee) 
Q. = pe WOS Wl 
Y Ro 


where o is the Stefan-Boltzman constant and Re represents 


the total resistance to radiation heat transfer between the 


plates as given by: 


Re = gee lb oe tee 
PL wo wl 


ie this form Eo represented the emissivity of the heated 
Temsheet surface and E74 was GHewemMissiviry or the adiabatic 
Plexiglas surface. 

Appendix B contains a representative heat balance and 


sample calculations for one of the experimental runs. 


Eee REOULTS 

Utilizing the expressions described earlier, the experi- 
mental data for buoyant forces aligned with the main direction 
of the fluid flow was reduced in accordance with the sample 
calculation procedure shown in Appendix B. Results for the 
major variables involved are provided in tabular and graphical 
form on the following pages. Table I and II show computed 
values for the straight and the curved test sections respec- 
tively. Figure 8 contains a plot comparing the average 


Nusselt number versus the Reynolds number for both test 
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Beerrons. Representative uneertainty bands have been indi- 
cated for all plotted experimental data and a complete error 


analysis 1s provided in Appendix A. 
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TABLE I 


Sv eee neneoUlls FOR STRAIGHT TEST SECTION 


_Re n(W/m2°C) _ Na. 
989 en 5.321 
876 19 25 5.298 
8 26 21.79 5.198 
751 Dale 5.197 
689 One u.Qud 
632 20.95 4.897 
560 19.20 u.559 
uae 18.79 4.482 
438 17.75 u.234 
363 16.99 u.053 
313 16.24 3.87H 
263 14.30 ecw 
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abe. Lt 


Set osseous TOR CURVED TEST SECTION 


_Re _De A(W/mec) =u 
eo) L44 2 mes. 6.694 
865 126 27.74 Ono ky 
SIRS alee, 25 06 Soe 
74L 108 2 one 2 De oO.0 
67S aS 22 ao o det 
654 Se 22a03 Spo on 
66 oid) Zora! zo s 
279 84 Zeb] Seog 
9904 oul 2134 Sd | 
20 76 20a u.875 
492 TZ ZOR Zs 4.838 
u29 6 3 ZOmec0) 4.819 
4O4 aC 2043 1. 4.844 
So 7 94 Rees u.,515 
342 90 eH ue 36 
Sls WS ibaa gale ora 
274 LO 14.50 3.460 
244 oo tol 2 Seuns 0 
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In all cases the results indicated an increase in the 
heat transfer rates as the flow velocity, 1.e., Reynolds 
number, was increased. However, when comparing the average 
Nusselt numbers of the straight and curved sections two 
distinct regions became evident. At the highest values of 
Reynolds number, which corresponded to the highest volu- 
metric flow rates, it was observed that the values of the 
average Nusselt number for the curved test section were 
greater than those in the straight test section. At the low 
and intermediate values of Reynolds number, the data was 
less conclusive with regards to the comparative values of 
the average Nusselt number in each test section. This region, 
between Reynolds numbers of 200 and 400, corresponded to 
values of the Dean number between 30 and 60. As Dean discussed 
Mommas work Concerning flow instability in a curved channel 
[12], a small disturbance would not be amplified below the 
critical Dean value of 36. Above this critical value, his 
linearized prediction indicated the disturbance would increas- 
ingly amplify as flow velocity increased. 

Consequently, it was concluded that the Taylor-Gortler 
vortices are only initially developing in the low to inter- 
mediate Reynolds number region and have not yet reached suf- 
PECTeMmE DGCOpertion to have any significant effect on the 


mecaeEsanscter rpates in the curved section. 


US 





The actual presence and onset of the Taylor-Gortler 
vortices was independently verified through the observation 
Seethe liguid erystal isotherm distribution on the Temsheet 
surface in the curved test section. At the smallest values 
of Reynolds number, the liquid crystals were in essentially 
straight bands that were perpendicular to the flow direction 
(Figure 9). A similar pattern was observed in the straight 
pest section for the entire experimental range of Reynolds 
numbers. This pattern was thought to indicate that the flow 
was in essence two-dimensional. As the value of the Reynolds 
number was increased, liquid crystal bands orientated in the 
streamwise directions became visible at the trailing edge of 
the curved test section. At the onset of these streamwise 
bands the liquid crystal bands at the leading edge remained 
orientated perpendicular to the flow direction. As the 
Reynolds number was further increased, the bands which were 
orientated in the streamwise direction eventually became 
evident over the entire surface of the curved test section 
(Figure 10). These streamwise bands would indicate a transi- 
PEt roen tEme initial, stable, two-dimensional flow to a less 
stable, three-dimensional flow triggered by the continuing 
development of the Taylor-Gortler vortices. 

The simultaneous occurrence of the streamwise bands, 
explained by the presence of the longitudinal Taylor-Gortler 


vortices, and the increase in the average Nusselt number 


U6 
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Deemer tpt lOneon Liguzd Crystal 
Cover sander curved» este section 
for Re=800 (De=117) 
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M@ecitted the conclusion that the vortices did in fact en- 
hance heat transfer from the curved test section. 

Comparison of the average Nusselt numbers for the 
straight and curved test section indicated that values for 
feme ratio of Nu, ie Nu. could be separated into two cate- 
gories. For Reynolds numbers between 400 and 800, the heat 
transfer rate was increased by an average of eleven percent. 
For Reynolds numbers above 800, the neat transfer rate was 
increased by more than twenty-five percent. In an earlier 
investigation of the boundary layer along a concave wall, 
Kreith [4], reported that heat transfer rates along concave 
walls increased from twenty-five percent to sixty percent 
for Reynolds numbers, based on the hydraulic diameter, 
between 10° and Moe This would suggest that the two data 
points for Reynolds numbers greater than 800 in the curved 
test section could represent flow that has tripped and 
become turbulent, or at least entered a laminar-turbulent 
transition. One possible explanation could be that the 
Taylor-Gortler vortices have started to break down at this 
a eeamicave sciereroOre Initiated transition to the turbulent 
flow regime. 

Experimental data obtained for the curved test section 
with the buoyant forces opposed to the main direction of the 
Flow indicated a minimal difference in heat transfer rates 


from the curved channel. This small decrease was only 
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evident at very low volumetric flow rates and became negli- 


gible for Reynolds numbers greater than 500. 
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bee Uso FONS AND CONCLUSIONS 


Investigation of the measured wall and air flow tempera- 
tures indicated there was negligible difference between the 
unheated wall temperature and the fluid bulk temperature. 

In addition, energy balances such as the one detailed in 
Appendix B provided evidence that the radiation heat transfer 
from the unheated wall was minimal. Based on this information, 
it was assumed that the heat convected to the air flow was 
solely from the heated wall. 

As previously discussed, the high aspect ratio of the 
channel provided the basis from which the initial assumption 
was made to model the experimental apparatus configuration as 
infinite parallel plates. The experimental data was found to 
substantiate this assumption as the computed values of the 
average Nusselt number in the straight section asymptotically 
Ee eoachegerie staeconetical yalue of 53.383 ,for parallel plates 
with one wall at a constant heat flux and the opposite wall 
adiabatic. [28] 

With the formulation of the problem complete, the experi- 
mental data could now be compared with analytical solutions 
of the same problem; i.e., laminar flow between infinite 
parallel plates with one wall at constant heat flux and the 


Opposite wall adiabatic. In particular, the experimental 


ome 


results were compared with two analytic studies; one con- 
Bmeted by MeCuem, et al. [26], for infinite parallel plates, 
and the other conducted by WorsgeSchmidt [27] for tubes 

and annular passages. In the case of the annular passages 
the results were compared to the numerical predictions for 
an annulus with an infinite inner and outer radius of curva- 
ture, thereby approximating parallel plates. 

Before comparing the experimental data with the studies 
mentioned above, it would be advantageous to address the 
initial assumptions and review the nomenclature common to 
Nepnethose studies. First, the following restrictions are 
applied in the mathematical formulation of the problem: 

1. the velocity profile is fully established 


2. the fluid transport properties and density are 
assumed constant 


ee eine: CONGUCTIOnNers neglipable 
4, viscous energy dissipation is negligible 
5. eonditions are invariant with time (steady state) 


All of these restrictions have been accounted for in the 
design of the experimental configuration and/or verified in 
the experimental process. Next, a dimensionless axial length 


Goordinate was defined as follows: 


me x 
s Oh Pr Re 


92 


where x was the length of the straight test section (29.2 ecm), 
Biewas the Calculated hydraulic diameter (1.239 cm), and Pr 
Wee the Prandtl number (0.70). In these studies, the 

Reynolds number (Re) was calculated on the basis of the 
hydraulic diameter. The values of x* indicated that the 
development of the experiments would most likely occur in 

the thermal entrance region. Analysis of the temperature 
distribution on the heated wall versus the axial length 
verified the flow was indeed thermally developing. 

In comparing the experimental results to those reported 
by MeCuen, et al., it was necessary to integrate their 
expression of the local Nusselt number for small values of 
x; 
ile 


0.67095 xxi’ 3 


Nu = 
Coos 


ote 


= 27 (lee )x* 





wo 
in order to obtain the following expression for the average 


Nusselt number. 








= -3 2 Swi, 2/3 
Nine = la ies = 0267055 ‘1 i Ou xs ] 
ue) ot? WO 
oe 


In this expression on would represent the heat convected to 
the air from the inner Plexiglas wall and Qe would be the 
heat convected from the outer Temsheet wall. The average 


(loco leummummerswas then caleulated in accordance with the 


ae 





experimental model where ee Vosme dud mero Zero and coe was 
synonymous with oe 

The experimental data obtained for the straight section 
was then recalculated as a function of the hydraulic diameter 
and both the experimental and analytical results were plotted 
mim deure -1. Note that the analytic expression for the 
average Nusselt number was only plotted for values of the 
Reynolds number that corresponded to x* less than 0.042 as 
this represented the maximum value for which the analytic 
solutions would be valid. Neglecting those experimental data 
points that plotted outside this valid range, there was good 
agreement between the experimental and the analytic results. 

A similar process was followed to compare the experimental 
results with those reported by Worsge-Schmidt. Again the 
analytic values for the local Nusselt number were integrated 
to obtain analytic values for the average Nusselt number. 

In this case the experimental and analytic values of average 
Nusselt number were plotted versus the non-dimensional axial 
length coordinate x*. The results are shown in Figure 12. 
Similarly, there was good agreement between the experimental 
and the analytic results. 

It is noted that in both cases where the analytic and 
experimental results were compared, the experimental data 
plotted above the analytic predictions. This difference can 


be explained, at least in part, by the limitations that are 
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inherent in any experimental configuration. In this case, 
the presence of side wall effects and the inability to 
totally account for all the heat transfer processes and/or 
losses combine to produce a slightly inflated vaiue of the 
heat transfer rate. Consequently, the average Nusselt number 
meealso slightly inflated. 

Irrespective of the slight differences locally, the over- 
all trends were conclusive. The heat transfer rates were 
shown to increase with increasing values of Reynold number. 
The ratio of ae ie ae indicated that there was a measureable 
increase in the heat transfer rate from the curved section 
once the Taylor-Gortler vortices had developed to sufficient 
strength. Although much less conclusive from this experi- 
mental data, it appeared the the Taylor-Gortler vortices 
were additionally important in the transition from laminar 
Eemeurbulent flow. 

The experimental data from the curved test section when 
the buoyant forces were opposed to the main flow direction 
was essentially inconclusive. While the trend was qualita- 
tively what had been expected, that is, the average Nusselt 
number for the buoyant forces aligned with the flow was 
greater than that for the bueyant forces opposed to the flow; 
the difference was insignificant in light of the experimental 
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Vi ese COMMENDATIONS 


While the amount of work remaining in the experimental 
study of heat transfer in a curved channel has been reduced 
for laminar flows, a large amount still remains to be 
accomplished. Additional experiments should be conducted 
to refine the correlation for the average Nusselt number as 
Beene plonsos the Dean number. Further work should certainly 
be considered in the region where the laminar flow begins its 
transition to turbulent flow. Detailed study in this area 
could provide important data to determine the manner of 
relationship that exists between Taylor-Gortler vortices and 
laminar-turbulent transition. With this work it will be 
necessary to design and construct a larger experimental 
configuration capable of higher volumetric flow rates and 
increased AT through the test sections. 

Finally it would be advantageous to do more experimental 
work in the general area of liquid crystal thermography, and 
more specifically, the photography of these crystals and their 


aster i bution ; 
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The uncertainties for the major variables in the experi- 
ments were calculated in accordance with the method described 
Py oe Kime vend —. Me€lintech [31]. The estimates of the 
uncertainty in the measured quantities were made conserva- 
tively so that there was considerable confidence in the cal- 
culated uncertainties. As an example, the calculation of the 


uncertainty for the Reynolds number is given below. The 


Reynolds number was defined by the equation: 





The uncertainty in the mass flow rate was determined as 


no lows-: 
dm _ dp.2 dR.2 dtT.2 dQ 
meee RC? te 


The uncertainty in the pressure, the gas constant, the tem- 
perature, the reading of the volumetric flow rate, the 


height of the channel, the dynamic viscosity and the cross 


a 





sectional area of the channel were obtained from estimates as 








follows: 
S- = 0.00048, SS = 0.00002, SF = 0.00056, = = 0.02000, 
Gite es dA 
os = 0.02003, —S—= = 0.00223, and © = 0.02054. 
Mair C 


The uncertainty in the mass flow rate was: 


oe y (0.00048) 7+(0.00002) 7+(0.00056) -+(0.02000)*=0.02001 


and the uncertainty in the Reynolds number could then be 


calculated as: 


SRE = ¥ (0.02001) “+(0.02003) “+(0.00223) “+(0.02054) “=0.03505 


The values of the uncertainties for all variables are 


itsted below: 


Quantity Uncertainty 
- 0.02054 
Apr, 0.00266 
cla 0.00415 

jaa 
De 0.04042 
Ie 0.02149 
; 0.00038 

_aif 
Nu 0.02938 
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700.5383 
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APPENDIA 8B 


SAMPLE CALCULATIONS 


Figure 13 provides a representation of the control volume 
used in the energy balance on the straight test section. The 
sample calculations that follow indicate the procedure used 
to obtain the heat transfer components required for the 
energy balance and eventually the values of the average 
Nusselt number as a function of the Reynolds number for the 
air flow. A similar energy balance and procedure was used 


mom the curved test section. 


ENERGY BALANCE 





Figure 13. Energy Balance in the Straight Test Section 


OZ 





SAMPLE CALCULATIONS 


A. DATA 
Taoom = eg 0°C 
Q = 45% 
Vor ieee ole 
ee =) SOROS HV 
Ror S F220 24 OR 
Tey —emipogzem) G=23,07 2C 
Trysq ieee Nie = 9 39.54 °C 
Try99 =-5 leo mV@p= 9 S645 OC 
Twol an) loos my 
Two2 = 1.843 mV 
Two3 =e eo my 
Two4 = 2.014 mV 
iWOm=—00,47 °C 
Twoo = 2.144 mV 
Two6 =e my 
Two7 = 2.274 mV 
Two8g e203 Se ITY, 
Twil oe 7 Seam 
Twi2 = 1.083 mV 
Tw13 =e. o Oem y 
Twi4 Seer i ol my 
Twid = 1.183 mV mace Oe °° 
Twi6 =e 12 mV 
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Twi7 =e, 2 Loma 


Twi8 = vile? hon 
eee Seco b= eens Oc 
ORG = 22.53 °C = 74.35 °F 
A = 0.0742 m2 
PL 
A = 0.0016 m° 
Cc 
A Se aeccr mn 
= oO 
—_ MO GeD amoC 
= Oo 
ep e 1.0057 KJ/Kg °C 
Magn = «2-983 x 107° Kg/m-sec 
_ -2 ; 
Kong eee ee sms 10 W/m °C 
AXrye = 0.00635 m 
e =O 0 
WO 
2 = 0.40 
Wl 
P eG 
WOrW1 
5 eee Gomee or Gym 


TEMPERATURE CALCULALIONS 


1. Temperature Difference in Insulation (AT Ig? 


= = = = = o 
AT Ns Tins1 Tons2 39.54-33.45 on Goa 


Tee “hempe ra ture (T,) 
In * tout 23.07 + 34.09 


i-CD 5 


3. Mean Temperature Difference (AT) 


he =. Peo ts e7 5556) = 21.89 °C 
wo B 


Boge 6 


64 








eo. POWER CALCULATIONS 


1. Power Supplied (Qp) 


a Vor “H _ (2.583)(30.008) _ eee 
a iene et )CC~=‘“ CS 


2. Heat Lost Through Outer Plate (Q)4? 


; AT INS 6.09 


Rho > = (00635) (CH. 18x10=2)(0.0742)) 
Lo AXrno/ Krys Apt Otic (C4. 18x10 020742 


Ze 


3. Heat Radiated (Q.) 


HR aGesaelon Resistance (Ro) 





Ll-€ l-e . 
ae WO 4 ale n Wi 
Ro Apt Ewo Spi Fwo-wi “PL Swi 
- < [= i = Sige fst 
PL wo WL PL 


b. Heat Radiated (Qr) 


ie = ——S ae 
R 


aK 


Lb 


(5.669x107°)(0.0742)(323.63'-301.78") 


ee 9 


3.843 W 
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4, Heat Convected to Air (Qo3y? 


a. Corrected Volumetric Flow Rate (q) 


Oo 
Teypp (OR? 


ee ee Ry 
STANDARD 


Peron) 534.04 
~ 60 a 529.69 


p00 24 pee 


DeeDeomsuty (pe) 


ae 1.013x10° 
RTpy,;, | C2E7Y(23. 534273. 16 
3 
= 1.189 Kg/m 


c. Mass Flow Rate (m) 
m= qp = (0.0024)(1.189) = 0.0028 Kg/sec 


Gee teat Convected to Air Ce 


4 
ut 


a m Cpair (Toye Try? 


CONC 02s) (i Uber est. 09-23 .07) 


it 


0.03103 KJ/sec 


Sepa. 


66 





AVERAGE HEAT TRANSFER COEFFICIENT (h) 


ae ace 2 ones 
Apr Mis oe Gomer 42) 21.99) 


19.11 w/M2°c 


AVERAGE NUSSELT NUMBER (Nu) 


Pee C19 11) (0100635) 8 
NU = sir = 0.02662 = 4.0959 


REYNOLDS NUMBER (Re) 


_ m d Pero O28 COR MeGs S ) 
NO = 4 ie = eee ee 
Hsin “ec (CGioee chm 9) Como 6) 
= 560.39 
= 560 


67 





re. 


Eee 


Pio imOr SRErFERENCES 


Pioneers Stability Of a Viscous Liquid Contained 
Between Two Rotating Cylinders", Philosophical Trans- 
Pommonsmen tie enhovale soclety of London, series A, V.223, 
Dpmegeg—-343, 1923. 


National Advisory Committee for Aeronautics, Technical 
Memorandum 1375, On the Three Dimensional Instability of 


Laminar Boundary Layers on Concave Walls, by H. Gortler, 
1942. 


Smith, A.M.O., "On the Growth of Taylor-Gortler Vortices 
Along Highly Concave Walls", Quarterly of Applied Mathe- 
matics, V.8, pp. 233-262, November 1950. 


Compas the Influence of Curvature on Heat Transfer 
tomuneamoressiple Friluids”, Trans. ASME, V.77, pp. 1247- 
2 Oe 95S. 


Mori, Y., and Uchida, Y., "Forced Convective Heat 
Transfer Between Horizontal Flat Plates", International 
Journal of Heat and Mass Transfer, V.9, pp. 803-817, 
1966. 


Tobak, M., "Hypothesis for the Origin of Cross-Hatching", 
Pievemommenail, V.d, No. 2, pp. 330-334, February 1970. 


Men piomeeemee Mopper, §.C., Blair, M.P., and Bailey, D.A:, 
"Fffect of Streamline Curvature on Film Cooling", 
leorenalwot Emgtneering tor Power, Irans. ASME, V.99, 
bemees SeNO. 15 pp. /7=382, January 1977. 


Nicolas, J. and LeMeur, A., "Curvature Effects ona 
itimweinerslade Cooling Film", ASME Paper No. /4-GT-156. 


Folayan, C.0. and Whitelaw, J.H., "The Effectiveness of 
Two-Dimensional Film-Cooling Over Curved Surfaces", 
ASME Paper No. 76-HT-31. 


Lord Raleigh, "On the Dynamics of Revolving Fluids", 
Proceedings of the Royal Society of London, series A, 
(oer elec oue 195. Reprints in Scientific Papers, 
V.6, pp. 447-453. 


ia emeeetl, Distribution ef Velocity and Temperature 

Between Concentric Rotating Cylinders”, Proceedings of 

the Royal Society of London, Series A, V.151, pp. 494- 
Zens 


68 





i. 


nS 


14. 


Le 


mG . 


me. 


ro. 


Eee 


Zi: 


Zee 


Dee 


Deans WeR., ‘Fluid Motion in A Curved Channel", Pro- 


eeeeeignmoueenec hoyal Society @f London, Series A, 
ar ie ape O2—420% 1928. 





LCV OnethesotablLlity of Viscous Flow in a 
Curved Channel", Proceedings of the Royal Society of 


London, Series A, V.244, pp. 186-198, 1958. 


Semeremranc. H., Boundary Layer Theory, /th ed., pp. 
529-536, McGraw-Hill, 1979. 


Mewmener si) ., Flentie, D.L., and McKee, R.J., "An 
Experimental Study of the Secondary Flow in a Curved 
Rectangular Channel", Journal of Fluids Engineering, 
VeweO epee os 2-360, March 1980. 


Viner ost. DUrao, D.F.G., and Crane, R.I., "Measure- 
ietEitebimecorhtler Yortices”, Journal of Fluids En- 
gineering, V.101, pp. 517-520, December 1979. 


Aihara, Y., "Nonlinear Analysis of Gortler Vortices", 
Themen vetesmonnrlumds, V.19, pp. 1655-1660, November 
eee 


oe eed Research Laboratories Report ARL 65-68, A 

SUpieeeoemewerodch tO the Enfluence of Gortler- -Type Vorti- 
Ceo Ouvne heat iransfer from a Wall, by Leif N. Persen, 
io 


McCormack, P.D., Welker, H., and Kelleher, M.D., "Taylor- 
Gortler Vortices and Their Effect on Heat Transfer", 
lepmenetmemmieatt Iranster,mimmg2, pp. g@1-112, February. 

ILS AT 


Kahawita, R. and Meroney, R., "The Influence of Heating 
on the Stability of Laminar Boundary Layers Along Con- 
cave Curved Walls", Journal of Applied Mechanics, 
WECM ll-l/. Marchilo7/). 0.  . 


Akiyama, M., Hwang, G.J., and Cheng, K.C., "Experiments 
on the Onset of Longitudinal Vortices in Laminar Forced 
Convection Between Horizontal Plates", Journal of Heat 

Transfer, V.93, pp. 335-341, November 1971. 


Ciegeweieec. did Akiyama, Mo, “Laminar Forced Convection 
Heat Transfer in Curved Rectangular Channels", Interna- 
tional Journal of Heat and Mass Transfer, V.13, pp. 471- 
ies Oe 7/0) 


69 





aS 3 


att. 


23. 


HO. 


27. 


iG 


Dai 


SHOR 


See 


Yee, G. and Humphrey, J.A.C., "Developing Laminar Flow 
and Heat Transfer in Strongly Curved Ducts of Rectan- 
gular Cross Section", ASME Paper No. 79-WA/HT-1L5. 


owe a oehi dca ana Ukon, D., "Forced Convective 
Heat Transfer in a Curved Channel with a Square Cross 
Section", International Journal of Heat and Mass Trans- 
Soe PD eo SCD 19/1. 


Durao, M. do Carmo, Investigation of Heat Transfer in 
Straight and Curved Rectangular Ducts Using Liquid 
Crystal Thermography, Eng. Thesis, Naval Postgraduate 


Pei oomnn Onteneygy Calitornia, 1977 . 


Department of Mechanical Engineering, Stanford Univer- 
sity, Report No. AHT-3, Heat Transfer with Laminar and 
Turbulent Flow Between Parallel Planes With Constant 
aicmemetaple Wall Temperature and Heat Flux, by P.A. 
NeCuewewene Kayse and WC. Reymolds, 12 April 1962. 


Wors¢ge-Schmidt, P.M., "Heat Transfer in the Thermal 
Entrance Region of Circular Tubes and Annular Passages 
with Fully Developed Laminar Flow", Internation Journal 
of Heat and Mass Transfer, V.10, pp. 541-551, 1967. 


Shah, R.K. and London, A.L., Laminar Flow Forced Con- 
vection in Ducts, Supplement 1, pp. 305-312, Academic 
ppece, 097s: 


Department of Mechanical Engineering, University of 
California, Berkeley, Numerical Computation of Buayancy- 
Induced Recirculation in Curved Duct Laminar Flow, by 
cmt and J.A.€s Humphreys March 1980. 


Poopecemince. sb ield, Kw, and Meyer, J.F., “Liquid 
Crystal Thermography and It's Application to the 
Study of Convective Heat Transfer", Journal of Heat 
imectingter. V¥.97. pp. 442-450, AWweust 1975. 


Memowoed. ang McClintock, FAs, "Describing Uncer- 
tainties in Single-Sample Experiments", Mechanical 


Emiimeeming, V./5, pp. 3-8, Janvary 1953. 


70 








Pee eVistTRIBUTION LIST 


Defense lechnieal Entormation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 


Associate Professor M. Kelleher, Code 69Kk 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93940 


Lieutenant James C. Ballard, III, USN 
United States Department of Energy 
Naval Reactors Representative's Office 
General Delivery 

Naval Base Branch Post Office 
Chagileeton-ssouth Carolina 29408 


Department of Mechanical Engineering, Code 69 


Naval Postgraduate School 
Monterey, California 93940 


Gil 


No. Coples 


























i iti 


